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We study a translation invariant spin model in a three-dimensional regular lattice, called the 
cubic code model, in the presence of arbitrary extensive perturbations. Below a critical perturbation 
strength, we show that most states with finite energy are localized; the overwhelming majority of 
such states have energy concentrated around a hnite number of defects, and remain so for a time that 
is near-exponential in the distance between the defects. This phenomenon is due to an emergent 
superselection rule and does not require any disorder. An extensive number of local integrals of 
motion for these finite energy sectors are identified as well. Our analysis extends more generally to 
systems with immobile topological excitations. 


I. INTRODUCTION 


Recently there has been a significant amount of interest 
on the mechanism behind how isolated quantum systems 
thermalize. While the eigenstate thermalization hypoth¬ 
esis asserts that quantum many-body systems thermal¬ 
ize at the level of the eigenstates,^’^ there are well-known 
counterexamples to this proposal as well. In particular, 
a seminal work of Anderson shows that disorder can in¬ 
hibit thermalization for noninteracting systems, as par¬ 
ticles are localized indefinitely at some fixed location.^ 
Recent studies indicate that the effect of localization can 
persist in the presence of interaction, even at finite energy 
densities.These results suggest that disorder can cause 
the system to “freeze” in time, and hinders equilibration. 
There are also some translation invariant systems which 
were proposed to be localized under a generic disordered 
initial conditions.However, currently there does not 
seem to be a decisive consensus as to whether such effects 
persist in the thermodynamic limit; see Ref. 14. 

In this paper, we point out that interacting quantum 
many-body systems can be localized by a completely dif¬ 
ferent mechanism. This is due to an emergent superselec¬ 
tion rule, and does not rely on disorder. We consider an 
exactly solvable and translation-invariant lattice model 
of spins, called the cubic code model^®, and find a set 
of states with manifestly localized spatial energy profile. 
These states remain invariant under Hamiltonian evolu¬ 
tion, and constitute the overwhelming majority of the 
states at finite energy. Our main result here is that all 
these properties remain stable against arbitrary pertur¬ 
bations, so long as the perturbation strength is smaller 
than some finite critical value. Hence we conclude that 
this phenomenon is a robust property of the phase. 

To be more precise, we consider a weak arbitrary per¬ 
turbation, denoted as T, to the cubic code model Hq that 
shall be defined later, and construct a set of states lify) 
such that their energy is concentrated around a finite 
number of points with respect to the perturbed Hamil¬ 
tonian Hq -P Y. These states are quasi-eigenstates in the 


sense that 
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where a,c,r] > 0 are constants and L is the smallest 
distance between the defects. Furthermore, the fraction 
of these states within the low-energy subspace is shown 
to approach I in the thermodynamic limit. 

While we only consider a system with finite number of 
defects, we emphasize that our system is fundamentally 
different from noninteracting fermionic systems. The 
transport of the low-energy excitations are hindered by 
a novel form of emergent superselection rule, which does 
not arise in free systems. Interestingly, localization orig¬ 
inates from strong interaction within our setup, contrary 
to what is observed in the context of disordered systems. 

It should be noted that the localized states satisfy¬ 
ing Eq. (1) are unlikely to exist in a system with non¬ 
trivial dispersion relation. Consider a one-dimensional 
fermionic chain with a Hamiltonian H = J dk e{k)clck as 
an example. Here, Ck is the annihilation operator of mo¬ 
mentum k. If the dispersion relation is e{k) = — Acosfc, 
then a state representing a single localized excitation, 
such a,s \ijj) = J -^^Ck l^)> spread out quickly in 
time; 


{iP\e IV”) = Jo(At) = 1 - ^ -h • • • , 

where Jo(t) is the Bessel function. 

In the rest of the paper, we shall introduce our model 
and explain how Eq. (1) is derived. We first begin by 
explaining the central concept, the locally gapped state. 
Roughly speaking, the notion of locally gapped state for¬ 
malizes the intuition that two eigenstates do not mix with 
each other with respect to a perturbation if they are ei¬ 
ther (i) separated from each other with a large energy 
gap or (ii) the corresponding off-diagonal matrix element 
of the perturbation is sufficiently small. We then pro¬ 
ceed by identifying all the low-energy states of our model 
which are locally gapped. A simple estimate shows that 
locally gapped states form the overwhelming majority of 
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the states within the low-energy subspace. Lastly, we 
rigorously prove that the aforementioned properties re¬ 
main stable in the presence of weak perturbations, and 
find an extensive set of approximate local integrals of mo¬ 
tion. Throughout the paper, we shall adopt the following 
technical terminology: A function f{x) rapidly decays if 
it scales like x°‘ exp(—for large x, where c, /? > 0 and 
a are constants. 
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II. LOCALLY GAPPED STATES 

With respect to a Hamiltonian H, a state \ijj) is said 
to be locally gapped with a diameter d and an energy A, 
if for any orthogonal state \(j)), one of the following holds: 

. MH\cj,)-{^\H\^)\>A>0 

• {(j)\ O I'(/')= 0 whenever O is supported on a ball of 
diameter d. 

This means that \ip) is separated by an energy gap of 
A from all the other states that are reachable by local 
operators from |'0). Unless specified otherwise, we will 
be interested in locally gapped states with an extensive 
diameter and finite energy gap. Let us mention some 
examples. A non-degenerate ground state of a Hamilto¬ 
nian with an energy gap is locally gapped. A less trivial 
example is the ground state sector of topologically or¬ 
dered systems. The ground states can be only connected 
to each other by an operator whose support is extended 
across the system. 

For clarity, let us also discuss states which are not lo¬ 
cally gapped. Consider a Hamiltonian H = ^ - erf on a 
chain where cr* is the Pauli-z matrix. It is gapped with 
a unique ground state in which all the spins are point¬ 
ing “down.” Its first excited state is a spin configuration 
with a single spin, say the jth one, pointing up. This ex¬ 
cited state is not locally gapped because there is another 
excited state for which only the (j -I- l)th spin is pointing 
up. These two are related to each other by a local opera¬ 
tor, namely crfcrf+i. These states are not locally gapped 
with diameter d > 1. 

The above conclusion exemplifies a general observa¬ 
tion: If the hopping of a quasi-particle is realized by a lo¬ 
cal operator, states with localized excitations are not lo¬ 
cally gapped. This intuition is true for many translation- 
invariant systems, and it is thus tempting to conjec¬ 
ture that the only way to have a locally gapped excited 
state is through some strong disorder. However, we shall 
show that there exists a translation-invariant system with 
many locally gapped excited states. 

One might want to consider an approximate and more 
general version of locally gapped states by allowing the 
transition amplitude (0| O to be approximately equal 
to zero for local operators O. We shall encounter this 
generalization later, but for now it will only make the 
notation more cumbersome. 


FIG. 1. Cubic code model. X and Z represent Pauli matrices 
and a^, respectively. A term in the Hamiltonian is the 
product of 8 Pauli matrices arranged as in the diagram. 


III. LOCALLY GAPPED STATES IN THE PURE 
CUBIC CODE HAMILTONIAN 


The cubic code model is defined on a simple cubic lat¬ 
tice with two qubits per site, where the Hamiltonian is a 
translation-invariant sum of the two types of interaction 
terms^®. 


Ho = -jY,Gf 
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See Fig. 1 for the arrangement of the Pauli marices in the 
terms. Every term in the Hamiltonian commutes with 
each other, and as such, the energy spectrum is discrete. 
For the ground state, both Gf and Gf take eigenvalues 
of -hi. The ground state subspace is degenerate and is 
topologically ordered in the sense that any local operator 
has the same expectation value for all the ground states. 
According to our definition in the previous Section, the 
ground state is locally gapped by an energy A = 2J 
and diameter L — 1, where L is the linear system size. 
The excited states can be described by defects, which are 
the violated local terms, e.g., Gf = — 1 or Gf = — 1. 
We refer to each types of the violated terms as A-type 
and 2-type. Since A-type term and 2-type term are 
lattice inversions of each other, it suffices to analyze just 
one of the two types. Not every configuration of defects 
is physically allowed. Indeed, a configuration of defects 
can be realized as an excited state if and only if it can be 
created by applying a tensor product P of Pauli matrices 
(Pauli operators) to the ground state. We call such defect 
configurations to be valid. 

The most general energy eigenstate is a superposition 
of valid configurations with the same number of defects. 
Note that the configuration of defects does not uniquely 
determine the state. For example, the ground state has 
no defect, but is degenerate. Given the defect configu¬ 
ration, this is in fact the only residual degeneracy. We 
shall call this as the topological degeneracy. 

Below, we will make use of the following fact^®. Under 
periodic boundary conditions there is an infinite family 
of system sizes such that any defect configuration with 
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an even number of X-type defects and an even number 
of Z-type defects is valid, and the topological degeneracy 
is 4. This is the family that we shall study in this paper. 
Concretely, it includes the lattices of linear size L = 2" + l 
for any integer n > 1. The situation is very similar to 
the toric code Hamiltonian; one cannot have an isolated 
anyon in a closed finite system, but any even number of 
anyons is allowed. 


A. No-strings rule 

The most important property, called no-strings rule^"^, 
is a formalization of the fact that isolated defects are 
immobile. More precisely, suppose \ip) is a state with a 
defect at site i and no other defects within distance d 
from i. If |'0') is another state with a defect at 7 ^ i and 
no other defect within distance d from j, then the no¬ 
strings rule asserts that ('(/''I O |^) =0 for any operator 
O supported on a ball of diameter d. Note that the no¬ 
strings rule is a nontrivial condition, in that i' can be 
even in the vicinity of i. 

The cubic code satisfies the no-strings rule, as we ex¬ 
plain below. Note that any operator can be written as 
a linear combination of tensor products of Pauli matri¬ 
ces (Pauli operators). This is possible simply because 
Pauli matrices generate the full operator algebra. This 
observation eases the problem; in order to argue that no 
operator can “transport” a defect, it suffices to show that 
no Pauli operators can transport the defect. If a Pauli 
operator O supported on a ball of diameter d is capable 
of “transporting” a defect within that ball, some product 
of O and its translations can transport the defect arbi¬ 
trarily far. The calculation of Ref. 15 implies that, if a 
defect can be transported arbitrarily far, it can be equally 
created from the ground state by acting some Pauli op¬ 
erator in its vicinity; see Fig. 2. This is impossible since 
a local Pauli operator can only create an even number of 
defects. We arrive at a contradiction, so we conclude that 
any Pauli operator, and hence any operator, is incapable 
of transporting defects. 



FIG. 2. No-strings rule. If some operator supported on a 
bounded region of the lattice can transport a defect, then 
there exists another operator that has exactly the same effect 
on the state and is supported on the vicinity of the initial and 
final defects. This is impossible since a single defect cannot 
be locally created or annihilated. Therefore, no such operator 
exists. See Ref. 15 for more detail. 


The no-strings rule implies the existence of low-energy 
excited states which are locally gapped. To see this, con¬ 
sider an excited state \ip) describing a configuration of 
defects that are separated from one another by a dis¬ 
tance d. Let e be the the total number of defects, so 
that the energy of the state is 2Je. If another state of 
the same energy, denoted as |(()), has a configuration of 
defects different from that of \tjj), there must be a de¬ 
fect at i that is present in \(j)) but not in By the 
construction of \il)), the no-strings rule implies that the 
matrix element ((/)| O |'0) is zero for any operator O whose 
support is contained in a ball of diameter d. More gen¬ 
erally, any state orthogonal to {ip) with the same energy 
can be labeled by two components. One represents differ¬ 
ent defect configurations with the same energy, and the 
other represents topological degeneracy. The transition 
from \tjj) to any orthogonal and topologically degenerate 
state requires an operator whose support is comparable 
to the system size. By the no-strings rule, the transition 
from 1 ^) to any other state with a different defect con¬ 
figuration requires an operator whose support is at least 
d. Other states have energy that is different from 2Je by 
at least 2 J. Therefore, the no-strings rule implies that a 
configuration of separated defects by distance d is locally 
gapped with diameter d and energy 2 J. 

Remark that the existence of the locally gapped ex¬ 
cited states is not due to any disorder; it follows from the 
no-strings rule, which is a nontrivial property of the cu¬ 
bic code model. A correct interpretation should be that 
defects at different locations represent distinct superse¬ 
lection sectors. A similar yet weaker version of this phe¬ 
nomenon is observed in Wen’s plaquette model^*. There, 
a defect can be transported to next-nearest neighboring 
site by a local operator, but not to the nearest neigh¬ 
boring site. In other words, the superselection sector 
is changed under a unit translation, although it is not 
changed under two units of translations. In the cubic 
code model, there is no such finite periodicity, and conse¬ 
quently, there are infinitely many superselection sectors. 


B. Typicality of locally gapped states 

Now we show that almost all defect configurations are 
locally gapped. Without loss of generality, let us assume 
that m defects are all A-type, as the analysis for the Z- 
type defects are exactly the same. Any configuration is 
valid as long as m is even, and there is no restriction on 
the positions of the defects. The total number of states 
with a fixed energy will be given by the number of distinct 
configurations, multiplied by the topological degeneracy 
4. This is a purely combinatorial problem in which one 
chooses TO defect locations from a lattice of V sites such 
that every defect is separated from the others by distance 
d. Let us imagine that we are placing one defect at a 
time, eventually to defects in total. The first defect can 
be anywhere, as there are V possibilities. The second 
one can be anywhere except near the first one; there are 
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V — V possibilities, where v is the volume of a ball of 
radius d. The third one can be placed anywhere among 
at least V — 2v positions. Continuing in a similar way, we 
see that the number of these sparse defect configurations 
is at least V{V — v){V — 2v) ■ ■ ■ {V — {m — l)v)/m\. The 
total number of all configurations of defects is (^). Thus 
the fraction of the locally gapped states with diameter d 
is at least 

V{y — v) ■ ■ ■ iy — {m — l)u) ^ 

The fraction approaches 1 algebraically in the system 
size, provided that d ~ L^~'^ where 0 < e < 1. We 
conclude that the majority of the defect configurations is 
locally gapped. 

There are excited states which are not locally gapped, 
but they form a vanishing fraction of the low-energy sub¬ 
space due to the previous analysis. These are the excited 
states with locally created clusters of defects. A local 
operator creates at least 4 defects from the ground state, 
and this cluster of 4 defects can hop to a nearby position 
by local operators. Thus, excited states with a topo¬ 
logically neutral (i.e., locally created) cluster of defects 
are not locally gapped, and may hybridize to become a 
momentum eigenstate upon perturbations. 


IV. LOCALLY GAPPED STATES UNDER 
PERTURBATIONS 


A. Band-preserving perturbations 

Consider a one-parameter family of Hamiltonians of 
the following form 

Hs = Ho + sY, se[ 0 ,l], 

where Hq is the original cubic code Hamiltonian (Eq. 2) 
and Y is some perturbation. We assume that the per¬ 
turbation consists of local terms with bounded strength, 
where strength of the perturbation refers to the supre- 
mum of the norms of the local terms. 

Let us invoke a powerful theorem which is due to 
Bravyi, Hastings, and Michalakis^®’^°. Setting the en¬ 
ergies of Hq to be En {n = 0 , 1 , 2 ,...) with the ground 
state energy being Eq = 0 , the theorem says that the 
spectrum of Els is contained in [E„(l — cs), A„(l -I- cs)] 
up to a small correction that vanishes in the thermo¬ 
dynamic limit. Here, the number c > 0 is a constant 
that is independent of the system size, and the strength 
of the perturbation is assumed to be smaller than some 
constant. At finite energies, this implies that there are 
energy bands that are separated from each other by some 
constant, say A. This spectral gap condition allows us 
to unambiguously define projectors onto the linear span 
of states within one of the separated bands. We shall 
refer to such projectors as the projector onto the band 
subspace. 

B. Evolution of the band subspace 


Locally gapped states are resilient to perturbations. 
Before we see why, let us recall a basic fact from per¬ 
turbation theory to develop some intuitions. Consider 
a nondegenerate eigenstate {ip) of some Hamiltonian H 
with an energy E. Under an infinitesimal perturbation 
Y, the eigenstate is perturbed as follows: 


The existence of the separated energy bands at finite 
energy implies that we can use a powerful machinery 
of quasi-adiabatic continuation^^’^^. This allows us to 
systematically track the effect of the perturbation on 
these subspaces. It says that, whenever a band subspace 
around energy is separated from the rest of the spec- 
trum by a constant gap, the projector Pi onto the band 
subspace obeys the following differential equation: 

a.pw = *[A(s),pW], (4) 


where |'0j) are the energy eigenstates of H with energy 
Ej. The correction term does not blow up if the state 
is either separated from the rest of the spectrum by a 
constant energy gap or the off-diagonal matrix elements 
are sufficiently small. These are exactly the two criteria 
that appear in the definition of locally gapped states. If 
the first criterion is satisfied, the state evolves under an 
adiabatic evolution, and we expect the wave function to 
change smoothly. If the second criterion is satisfied, the 
perturbation acts almost diagonally with respect to the 
original energy eigenstates, so we expect the effect of the 
perturbation to the eigenstates to be negligible. These 
observations suggest that if the energy profile of \ip) is lo¬ 
calized in the sense that there are a few positions around 
which energy density is nonzero, then the perturbed state 
also has a similar energy profile. This rough intuition will 
be made more precise in this Section. 


for some hermitian operator K{s). Formally, the solution 
to this equation is the following: 

py = UsPt^u! (5) 

Us = T exp ^ J K (m) drt^ (6) 

where T is the “time”-ordering operation. The readers 
may wonder about the definition of K{s) at this point. 
While its rather esoteric definition shall be given in the 
Appendix, its essential property can be summarized suc¬ 
cinctly as follows. 

The important point is that the unitary Us preserves 
locality. Namely, for any operator O supported on a re¬ 
gion M, we have the following decomposition: 

yous = Y.op, iio;ii </(r), 


( 7 ) 
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where is some operator supported on a distance-r 
neighborhood of M and /(r) is some function which de¬ 
cays rapidly in r. 

Note that Us governs the evolution of any energy band 
which is separated from the rest of the spectrum by a 
constant gap. As we have already discussed, the energy 
bands of the perturbed cubic code at finite energy sat¬ 
isfies this condition, provided that the strength of the 
perturbation is sufficiently small. 

C. Approximately locally gapped states 

By making use of these results, we can explicitly 
construct a large set of (approximately) locally gapped 
states. Such states are defined as follows: 

\ips) = Us IV'o), (8) 

where IV'o) is any state describing a sparse defect config¬ 
uration of Hq, the unperturbed cubic code Hamiltonian. 
Let us assume that (i) the perturbation strength is suf¬ 
ficiently small so that I'i/'o) belongs to some separated 
band, say and (ii) IV'o) is locally gapped with a 

diameter d. 

These states are approximately locally gapped in the 
following sense. For the states that are not in the band, 
there is an energy gap A. For the states in the band, 
we need to consider the transition amplitude. Suppose 
there is an operator O supported on a ball of diameter 
d/2 such that (0s| O IV's) = 5 where {(j)s\'4’s) = 0. This is 
to say that ((/)o| U/OUs \ipo) = ^ where \(j)o} = U/ \(j)s) is 
a state that belongs to Using the decomposition 

in Eq. (7), the operator U/OUs is quasi-local and hence 
can be approximated by an operator O' supported on 
the ball of diameter d such that WU/OUs — 0'\\ decays 
rapidly in d. Since |^o) is locally gapped with radius d, 
(0o|O'|'0o) = 0. It follows that 5 decays rapidly with 
d. This establishes that |'0s) is approximately locally 
gapped with a diameter of d/ 2 , up to an “error” that 
decays rapidly with d. 


V. PROPERTIES OF LOCALLY GAPPED 
STATES 

In the previous section, we have identified the locally 
gapped states of the perturbed Hamiltonian(iJs). In this 
section, we study their static and dynamical properties; 
we show that they have spatially localized energy profile, 
obey an area law of entanglement entropy, behave almost 
like energy eigenstates of Ug, and admit an extensive 
number of approximate local integrals of motion. Since 
the ratio of the number of localized states relative to 
the dimension of the low-energy subspace is very close 
to 1 from our counting argument, we conclude that a 
typical state from the low-energy subspace of arbitrarily 
perturbed Hamiltonian consists of “frozen” defects. 


A. Static properties 

Here, we show that the spatial energy profile of |'0s) 
is localized. We actually make a stronger claim: that 
the local reduced density matrix of |'0s) at a region M is 
close to that of the ground state of Hg if M is far from 
the locations of the defects for \tpo). 

In order to see this, recall that the local reduced den¬ 
sity matrix is completely determined by the expectation 
values of local observables. 

{^s\0\i;s) = ii^olU/OUsli^o). 

Since Us preserves locality, the “dressed” operator 
U/OUs can be well-approximated by some local opera¬ 
tor, say O'. If this local approximation acts away from 
the defects, its expectation value with respect to ji/o) re¬ 
produces the ground state expectation value of the unper¬ 
turbed system. If the ground states of the unperturbed 
and perturbed Hamiltonians are |Ho) and |Hs), respec¬ 
tively, the preceding discussion implies that (^s| O |V's) — 
(HolO'l^^o) = (Hg|C/gO'C/t |Og) ~ (Hg|0|Hg). There¬ 
fore, any observable acting far from the defects cannot 
distinguish lips) from jUs). The error term decays rapidly 
in the distance from O to the defect’s location. (Recall 
that a rapid decaying function scales like for 

some constant c, /3 > 0 and a for large x.) In comparison, 
it is worth noting that the difference between the local 
reduced density matrices of a momentum eigenstate and 
that of the ground state can only be algebraically small 
{x~"' for some 7 > 0) in the system size. 

We note in passing that the state \tps} obeys an area 
law of entanglement entropy. Namely, the entanglement 
entropy S{M, lips)) for a ball-like region M scales like 
the area of the boundary of M. This is not too sur¬ 
prising in light of a recent result^^; in fact, it has noth¬ 
ing to do with our construction of states with a local¬ 
ized energy profile. This is rather a consequence of the 
locality-preserving property of Us and the smallness of 
the low energy band subspace. Consider first the unper¬ 
turbed Hamiltonian Hq. Every state ji/o) with a defi¬ 
nite configuration of m defects, that is not necessarily 
sparse, has the same entanglement entropy for any re¬ 
gion as the ground state, since it is related to the ground 
state by {ipo) = P |Ho) for some tensor product unitary 
operator P. Moreover, since the ground state is de¬ 
fined by a commuting set of Pauli operators, the nonzero 
Schmidt coefficients with respect to the region M are all 
the same. There are precisely 2^° nonzero Schmidt coef¬ 
ficients, where Sq = S{M, \'ip q)) / \og2. We consider the 
Schmidt decomposition of an arbitrary superposition |^/o) 
of the states within the band Pg™^ with respect to the 
region M and its complement. Since there are at most 
C>(Vol M)™ states with different defect configurations 
in M, the Schmidt rank is at most ©((Vol 
The topological degeneracy does not enter to this up¬ 
per bound since topologically degenerate states have ex¬ 
actly the same reduced density matrix for a ball-like re¬ 
gion M . The entanglement entropy is therefore at most 
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S{M, I'i/'o)) + mlog Vol M + 0(1), which is an area law. 
The result of Ref. 23 states that for any state \tj}), the 
evolved state Ug IV') obeys an area law of entanglement 
entropy whenever |^) does. Since any state \ips) in the 
band is equal to Ug \ipo) for some |^o) in we 

have the area law for lips). 


B. Dynamical properties 


Now we are in a position to derive the main result, 
namely Eq. (1). The key insight is that the local Hamil¬ 
tonian has small matrix elements between the locally 
gapped state Itpg) and the other states. To formalize this, 
let us fix a basis of the band subspace consisting of 

|1) = IV's) = Ug IV’o), k) = Ug |g) (g > 1) 

where \q) are energy eigenstate of Hq with definite defect 
configurations. Here |i) is approximately locally gapped, 
but \q) may not be so. In the Appendix, we show that 

^MHg\-ips)\<g{d) (9) 

where g is a rapidly decaying function. Since the matrix 
element is so small, the time evolution can be “simulated” 
by neglecting those small matrix elements. Indeed, let Q 
be a hermitian matrix acting within PP , defined as 


{u\Q\v) 


0 if u = 1 ^ V or u ^ 1 = V, 

{u\Hg\v) otherwise. 


( 10 ) 


Then, Eq. (9) implies that 

\\PgiQ-Hg)Pg\\<gid), (11) 

where we used the fact that the operator norm is at most 
the sum of absolute values of the matrix elements. This 
in turn implies that 



f\t\ 

< dw ||e™'^(A - P)e-™-®|| < \t\g{d) 

Jo 


where A = PgHgPg and B = PgQPg, from which we 
conclude that 


|(V's| e \pjg) \ > 1 - \t\g{d). 

This is Eq. (1). 

The fact that these localized states remain almost in¬ 
variant under time evolution is closely related to the 
existence of the local integrals of motion. The defini¬ 
tion of local integrals of motion in our context is a bit 


more relaxed than those that are discussed in the con¬ 
text of many-body localization. Since we are studying 
a finite-energy subspace, the integrals of motion of our 
interest may not commute with the Hamiltonian on the 
entire Hilbert space. In fact, since our system can be 
translation-invariant, we should not expect such commu¬ 
tation relation in general. However, there are local op¬ 
erators that commute with our Hamiltonian within the 
localized subspace, e.g., the subspace spanned by the lo¬ 
cally gapped states of the perturbed Hamiltonian. 

More concretely, let us define Pioc be the projector onto 
the localized subspace. This refers to the linear span of 
Ug\q), where |g) describes a state with defects which are 
separated by a distance d. Define 

I, is) = Ughfu! (12) 

where is the local term of the Hamiltonian Hq. 
This operator is quasi-local due to the locality-preserving 
property, i.e., Eq. (7). 

Now let us see whether this operator commutes with 
the Hamiltonian within the localized subspace. 


\Pioc[Hg,I,{s)]Pioc\\ < ^ |(a|[i7„X,(s)]|6) 


a^b 




a,b 


<^2 {a\Hg\b) < 2(dimPioc)g(d) 

a^b 

< 2L^^g{d). (13) 


where L is the linear system size. The first inequality 
is the triangle inequality, the second one follows from 
ll^j(®)ll = Ij the third one is from Eq. (9). If we 
choose d ^ L^~'^ for a small e € (0,1), then the upper- 
bound decays rapidly in the system size, and Eq. (3) 
reads 


dimPi 

dim 

for some /3 > 0. Thus, Ij{s) is almost an integral of 
motion for almost all states of the band. Indeed, the ex¬ 
pectation value of Ij{s) with respect to locally gapped 
states is slowly varying. Let lip) be any state in the sub¬ 
space Ploc. 

(V^|P‘^=X,(s)e-‘«» IV^) 

= (V’|2(,-(s) IV') + f dw {'ip\P'"^H[Hg,Ij{s)]e~^^^‘ {ip) 
Jo 

= (V'|X,-(s) IV') ± f dwO (||Ploc[Ps,Xj(s)]Ploc||) 

Jo 

= (V'|X,(s)|V')±iO {L^9{d)) 

The expectation value of Bj{s) stays at the initial value 
up to a time that grows rapidly in the system size. 
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VI. DISCUSSION 

There is a considerable amount of work in the liter¬ 
ature regarding the fate of the disorder-driven localized 
system in the presence of interaction; see e.g. review ar¬ 
ticles Ref. 24 and 25. Here we have taken a different 
route, and identified a mechanism for localization which 
is driven by a strong interaction. By exploiting the fact 
that the ground state of our model is topologically or¬ 
dered, we rigorously showed that the localized energy 
profile of the majority of the low energy states remains 
unchanged under arbitrary perturbations. 

We remark that, though our analysis was focused on 
one particular model, a similar phenomenon is also ob¬ 
served in other models such as Chamon’s spin model^®’^^ 
and other “cubic code” models^®’^®’^®. In Chamon’s 
model, there exists a string operator that is capable 
of moving topologically nontrivial excitation. However, 
there exists an isolated excitation such that no opera¬ 
tor can move it without creating extra excitations. A 
sparse configuration of these isolated excitations is a lo¬ 
cally gapped state, and our analysis carries over. 

Our work strongly suggests that the prevalent dichoto¬ 
mous view on the role of disorder and interaction needs to 
be modified. That is, certain strong enough interactions 
can drive the system into an exotic topologically ordered 
phase, in which localization occurs at low energy. 

An outstanding open problem is whether a locally 
gapped state can exist at hnite energy density, and re¬ 
main stable against perturbations. At such high energy, 
the gap between the bands would collapse upon perturba¬ 
tion, and the machinery of quasi-adiabatic continuation 
in its present form cannot be applied. 
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Appendix A: Generator of quasi-adiabatic 
continuation 

The quasi-adiabatic continuation operator([/s) obeys 
the following differential equation: 

^Us=tK{s)Us, (Al) 

as 

where Uq = I. Here we explicitly write down the gen¬ 
erator, namely K(s). While there can be many different 
choices of the generator, the specific choice we make is 


due to Bachmann et al^^: 

/ OO nt 

w^{t) / (A2) 

-OO J 0 

Here w-y{t) is a normalized and rapidly decaying function 
whose Fourier transform is supported on [— 7 , 7 ]. Such 
a function is known to exist. An example is given in 
Ref. 30, which decays faster than exp(—cyt/log^( 7 t)) for 
some c > 0. If the energy gap between a band and the 
rest of the spectrum is bounded from below by 7 , K{s) 
generates the unitary flow of the band subspace along the 
path. The rapid decay of (t) implies that Ug preserves 
locality. In particular, it is known that Ug maps a local 
operator to a quasi-local operator, i.e., an operator that 
can be well-approximated by a local operator. 

Appendix B: Proof of Eq. (9) 

Consider the off-diagonal element of the Hamiltonian 
between two orthogonal states within the band subspace. 

{q\Hg\i) = Y,{q\u!h,Ug\i) 

i 

The term U^hiUg is quasi-local, and thus can be decom¬ 
posed as X]r>o r-! where h- ^ is supported on a ball of 
diameter r centered at i such that ||h'^|| < /(r) for a 
rapidly decaying function /. We need to consider two 
cases. 

Firstly, if |g) and |1) have the same defect configu¬ 
ration, i.e., they are topologically degenerate, then hi y. 
must have support across the system; r ^ L. The topo¬ 
logical degeneracy given a defect configuration is 4 for 
the cubic code under the aforementioned boundary con¬ 
ditions. Thus, the contribution from the topological de¬ 
generacy is at most 4/(L). 

Secondly, suppose |g) and |1) have different defect con¬ 
figurations. For each f, let j(i, q) be the furthest location 
from i at which the defect is present, only in one of the 
two states. Each off-diagonal terms, i.e., (gl h' ^ |1), van¬ 
ish unless h[ ^ creates or annihilates a defect at j = j{i, q). 
In order for this to happen, r has to be larger or equal 
to the distance between i and j. Therefore, if dist(f, j) is 
large, the off-diagonal terms decay rapidly. If dist(i, j) is 
small, a better bound can be obtained as follows. Since 
| 1 ) is locally gapped with a distance d, (g| h - ^ | 1 ) vanishes 
for r < d. These two cases cover all the possibilities, and 
we conclude that the nonzero contributions come from 
h' J. with r > R(i, q) = max(dist(f, j), d): 

|(g|C/th,C/,|l)| < ^ /(r)=:/i(i?(*,<z)). 

r>R{i,q) 

We now sum the right-hand side over all the states q and 
terms hi, to obtain Eq. (9). 

Suppose m is the number of defects in |1). Let 
us upper-bound the number of tuples (i, q) such that 


x3m+3_ 


dist(z, j(z, q)) < a; for a given x. The defects of Ig), as well 
as the site i, has to be contained in a 2a:-neighborhood of 
the defects of |1); otherwise the distance between i and 
j is larger than x. The number of such configurations is 


at most ((2a;)^)™ • m{2x)^ ^ 

i,q 

= 51 /i(^)+ 5] 

i,q:dist{i,j)<d i,q:dist{i,j)—x'>d 

<0(l)5]a:"/i(x), 

x>d 


where in the last line we made use of the fact that the 
topological degeneracy is finite. Accounting for the topo¬ 
logical degeneracy, we get 

5] |(g| |1)| < 4/(L) + 0(1) 5] x^^+^Mx). (Bl) 

q^l x'>d 

Both / and /i are rapidly decaying. Eq. (9) is proved. 
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